Exact Yrast Spectra of Cold Atoms on a Ring 
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We propose a methodology to construct excited states with a fixed angular momentum, namely, 
"yrast excited states" of finite-size one-dimensional bosonic systems with periodic boundary con- 
ditions. The excitation energies such as the first yrast excited energy are calculated through the 
system-size asymptotic expansion and expressed analytically by dressed energy. Interestingly, they 
are grouped into sets of almost degenerate energy levels. The low-lying excitation spectrum near 
the yrast state is consistent with the U(l) conformal field theories if the total angular momentum 
is given by an integral multiple of particle number; i.e., if the system is supercurrent. 
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Ultracold quantum gases have attracted great interest 
due to their potential application to the testing of the 
quantum fluctuations in many-body systems [ij. Non- 
trivial correlations in one-dimensional exactly solvable 
quantum many-body systems have been theoretically 
studied in the framework of the Bethe ansatz and con- 
formal field theories (CFT) with finite-size scaling 0, ■ 
While these systems are known to reveal large quantum 
fluctuations, there had been no ideal experimental test- 
ing ground until various control techniques in ultracold 
gas, such as changing the effective spatial dimension, sign 
and strength of interaction [if, became available. 

A periodic one-dimensional (ID) bosonic system, for 
example, Lieb-Liniger (LL) model [H, is simple mathe- 
matical model with various applications to the study of 
many-body phenomena in a wide range of interactions . 
The periodic geometry has been experimentally realized 
as a circular wave guide or elliptic trap 0, H[ ■ In periodic 
ID systems it is possible to study rotational properties by 
giving finite angular momentum to the system. Because 
the LL model has translation symmetry, it is a useful 
tool for studying the energy spectrum within the Hilbert 
subspace of a fixed angular momentum. This concept, 
originally called the yrast problem in nuclear physics [9j , 
is known to provide an insight into ensembles under ro- 
tation, the generation of topological defects 9 Ml3| . and 
the study of quantum Hall phase transitions [14J . 1 1 51 ] . In 
particular, two-dimensional yrast problem has been well 
studied in the context of vortex creation in a harmonic 
trap. The yrast state with the angular momentum of 
non-integral multiple of hN is regarded as an off-centered 
vortex state when it is viewed from the Gross-Pitaevskii 
mean-field picture, while that of ±HN is the centered 
single vortex state 10]. This approach would also give 
some important clues to understanding the fundamental 
properties of superfluidity 3- H| ■ 

In this Letter we construct low-lying excited states for 
the LL model with a fixed angular momentum. The 
yrast states in the LL model are known to correspond 
to the lowest-energy hole excitations, the set of which 



are called the type-II excitation branch in the thermo- 
dynamic limit jl9| . However, we note that the concept 
of the yrast problem is defined only in a finite-size sys- 
tem. We define, in this letter, the excitation energies 
from the yrast state as "excited yrast energies" . We an- 
alytically express the lowest-order yrast excitation ener- 
gies in terms of dressed energies through the asymptotic 
expansion with respect to the system size. We evaluate 
these excitation energies by obtaining the particle quasi- 
momenta as an asymptotic solution of the Bethe ansatz 
equations, and by defining the hole quasimomenta. Such 
excitations play a fundamental role in the spectral analy- 
sis related to scattering experiments [2(| when the system 
under probe has finite angular momentum. 

The LL Hamiltonian for N bosons interacting with 5- 
function potential in ID is given by 
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Here <?id is the ID coupling constant [21] . We consider a 
circular geometry with the radius R by imposing periodic 
boundary conditions. The coordinates Xj are given as 
Xj = R9j with 8j being the azimuthal angle < 9j < 2tt, 
and £ = 2ttR the circumference of the ring joining these 
coordinates. 

We first discuss the Bethe ansatz theory which pro- 
vides exact eigensolutions to the Hamiltonian (fTJ). The 
set of quasimomenta k/s for j S {1,2,..., N}, satisfies 
the Bethe ansatz equations (BAE) 



JY 



£k 3 = 2-nIj 



2 arctan 

otft 



(2) 



Here u = mgm/h 2 is the coupling constant with the 
dimension of wavenumber. Each eigensolution of the 
Hamiltonian ([1]) is specified by a set of quantum num- 
bers Ij, which are integers for odd N, and half- integers 
for even N. The ground-state solution, for instance, cor- 
responds to the set of integers 1^ = j — (N + l)/2 for 
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j G {1, 2, . . . , N}. It is convenient to define the maxi- 
mum and minimum quantum numbers (fermi points) for 
the ground state as 7* = ±(N — l)/2. By solving the 
B AE for a given set of Ij , we obtain the set of quasimo- 
menta kj . Then the eigenvalue of the Hamiltonian (TTJ) is 
obtained by summing up k 2 [f|. The total angular mo- 
mentum is given by L = RK XwLi kj, or ; according to 
the BAE, alternatively by 
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Hereafter we shall often denote L/h by L for simplicity. 
In the strong coupling Tonks-Girardeau (TG) limit u —> 
oo, the second term on the right-hand side of Eq. ((2]), 
approaches zero, and the solutions are, therefore, given 
by kjl = 2-irIj. 

In order to evaluate quasimomenta kj of an arbitrary 
eigensolution, we introduce a function 
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In terms of Zg(k), the set of BAE @ is simplified to 

Ze{k j ) = I j- for je {1,2,..., N}. (5) 

When N and I are finite but large, the summation of 
kj is approximated to an integral including a finite-size 
correction with the use of the Euler-Maclaurin expansion 
for large I. We investigate the approximate BAE solu- 
tions for the ground state and excited state of a pair of 
particle-hole excitation as follows. 

It is easy to show that the ground-state root density 
p(k) satisfies the integral equation 
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where K{k) = 2u/(k 2 + u 2 ). We have defined ±Qq by 
Z e (±Q Q ) = (I ± ±l/2)/£. They are related to the particle 
density by N/£ = J^q o p(k)dk. Approximating Zg(k) we 
define a ground-state counting function Z(k) as an inte- 
gral of the root density Z{k) = J_q o p(q)dq + Z(—Q ) 7 
where p{q) denotes a solution to Eq. ©. With Z(k) 
as a function of k, we have a mapping between kj and 
Z(kj) = Ij/£ via Eq. We evaluate kj's by the in- 
verse function of Z; i.e., kj = Z~ l (Ij/€). The difference 
between the (particle) quasimomenta, evaluated by the 
inverse of Z, and exact solutions, as direct solutions of the 
BAE, is of the order 0(l/£ 2 ). Here, we derive Z(k) from 
Zg(k) replacing the sum with the integral, which gives 
corrections of the order 0(l/£ 2 ). We thus denote the 
ground-state quasimomenta kj = Z~ 1 (7j 9 ' ) /£) by k[I^}. 



We derive a particle-hole excited state from the ground 
state making a hole at Ih and a particle at I p . We denote 
by kj the quasimomenta satisfying the BAE with Ij = 
I ( j 9) for 1 < j < N and j + j h {j h = I h + (JV + l)/2), and 
by k p that of I p . We define Zg(k) replacing kj's in Eq. 
(j4| with kj's and k p . We define the hole quasimomentum 
kh by Zi(kh) = Ih/£- In the ground state we define kh 
by Z({k h ) 
satisfies 



I h /£. The root density ^{k) = dZ £ (k)/dk 



Pf h (~k) 



1 

2^ 
f 

~~£ 



Q 
6{k 



dq 

- 2^ 



K{k 



l)p P i' h (l) 



1 



48tt£ 2 

-o(r 3 ) 



-k h )-S(k- 

' K'(k-Q+) 
, p p / 1 (Q + ) 



K'(k-Q-) 



pjl 



pV(Q-) 



Here Q ± are related to N and L bytfgl p p e ' h (k)dk = N 
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(4) and R£ kp p t ' h (k)dk = L, respectively. 



Shifting the Hamiltonian as T-Ll, n 
we define the dressed energy [3] by 
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where eo(fc) = (hk) 2 / (2m) — p — VLRhk. We determine p 
and VL so that the dressed energy vanishes at Q ± : 



e(Q+) = e(Q-) = 0. 
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In the asymptotic expansion with respect to £, we obtain 
by Eq. ([7]) the particle-hole excitation energy E^ h N as 
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The exact result generalizes both the type I and II exci- 
tation branches [19j, including the finite-size correction 
of the order 0{l/£). 

For a given angular momentum L with fixed particle 
number N we can construct the q-th eigenvalue E^ N 
systematically for q G N. We denote the particle-hole 
quantum numbers as 7^ for the q-th state. It is sufficient 
to consider only a domain < L < N due to the Galilean 
invariance. We often simplify E^ N as E^ ] . In Fig. Q] 

eigenvalues E^ are plotted, where the lowest line shows 
the yrast spectrum q = 1 as a function of L/N. 

Let us construct the yrast state q = 1 (type II branch 
found by Lieb in the thermodynamics limit 19]). It is 
derived from the ground state of L = by making a hole 
1 — 7 and a particle at Ip = I + + 1, as 



at 7« = 7 



shown in Fig. [5] (i). Here we remark that Tp 1 ' — I^ 1 ' = L. 



r(i) - 



We evaluate the quasimomenta by k 
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particles at {Ip 2 ^} 
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state. We denote the excitation energy associated with 
these "right" processes as AE^. We have 



= {Ip 1 ^ + 1, Ifr 1 '} and two holes at 
Ii.*' + 1} are excited relative to the yrast 



AE+ = e(fc (2) [7+ + 2]) - e(fc« [7+ + 1]) 

+e(fc (1) [I+ + 1-L])- e(fc (2) [7+ + 2 - L]). (12) 



FIG. 1: (Color online) Eigenvalues of LL Hamiltonian 

measured from the ground-state energy Eq for u = 1.0 (N = 
11) by solving the BAE. 



Several low- lying eigenvalues E L ^ measured from the 
yrast-state energy E^/ are plotted in Fig.[3J Each branch 
is specified by a particular set of particle-hole pattern, 
and is linear with respect to L/N. 

Let us calculate the first excitation energy from the 



In the TG limit, Eq. (HI is simplified to AE+/A = 27 
where A= (h 2 /2m)(2Tr/£) 2 . 

(b) Left excitation: move a particle at I~ to I~ — 1, and 

1. We thus have particles at 



move a hole at ij^ to ijp 



{if)} = {/- - 1, 1^} and holes at {7^} = {/", I, 
1}. The left excitation energy reads, 

A^= £ (fc( 2 )[7--l])-e(fcW[7-]) 

+e(k (1) [I + + 1 - L\) - e(k (2) [I + ~ L}), (13) 

which reduces to AE^/A = 2(N - L) in the TG limit. 

The first excitation energy AEl as a function of L is 
given by 



(2)! 



r(i) 



yrast state AEl = E^ 1 — E L ^\ By evaluating the quasi- 



momenta for q = 2 as fc) 
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7 = 0,1, 
2 < 7 < N/2, 
N/2 < L < N - 2, 
L = N-1,N 
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Here, the first term in Eq. (fT0| changes for the states of 
q = 1, 2 only by the order 0(l/£ 3 ) due to Eq. ©. 

The excitations from q = 1 —> 2 are constructed either 
of the following processes, as shown in Fig. [2J 
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where we regard the absence of a particle at I + + 1 in 
the ground state as a hole at 7^ 1 - ) for 7 = 0. Similarly for 
7 = N, we regard the absence of a particle at 7~ in the 

(a) Right excitation: move a particle at I ( p 1] to I p 1} + 1, ground state as a hole at 7^\ For 7 = 1 (7 = N - 1) 

it is obvious from Fig. [2] that the right (left) excitation 
is impossible because there is only one particle next to 
the hole at The AEl is symmetric with respect to 
7 = N/2, and minimal at 7 = 2 and 7 = TV — 2. 

Through the asymptotic expansion with respect to I, 
the left and right excitation energies AE^ are given in 
terms of the ground-state quasimomenta as follows. 



and move a hole at ii 1 ^ to I, 
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1. As a result, two 
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The difference between k^ and kj is of the order 0(1/ 1) 
for each j. It is also the case for fc/, and k%, and hence 
we have put the ground-state momentum k^ in Eq. (I15[) . 

Some branches with different particle-hole configura- 
tions are "quasi-degenerate" as shown in Fig. |3] For the 
nth lowest set of E L ^s, at most n — 1 energy levels are 
lying within the range of the order 0(1/^ 2 ). Here, the 
energies of higher excited states q > 3 are systematically 
calculated by giving particle-hole configurations similarly 
to the case q = 2. Moreover, for 7 close to the integral 
multiple of N, particle-hole excitations cost large energy 
due to tight constraints on configurations. 




0(1/?). (15) 
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FIG. 2: (Color online) The scheme for constructing the (i) 
yrast ground (q = 1) and (ii) first excited (q = 2) states 
for N = 5. The integers denote the quantum numbers Ij's. 
(i) A single hole is present at i£ . (ii) There are two types 
of excitations, where a particle and a hole move to the left 
(L — 0,1,3), and to the right (L = 2,4,5). The change in 
the energy associated with the left and right excitations is 
denoted as AE£ and AE~/, respectively. 
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FIG. 3: (Color online) Eigenvalues of LL Hamiltonian 
measured from the yrast state for u = 1.0 (N = 100) 
through solutions of the BAE, where E 



(q) 



for q 



2, 3, ■ ■ • , 7 correspond to blue, red, green, purple, orange, and 
light blue dots, respectively. The nth set of branches of E 
E^ has at most (n — l)-fold quasi-degeneracy. 
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In order to numerically implement the above method, 
we calculate p(k) and e(k) solving the integral equa- 
tions dU) and dl]), respectively, setting Q ± — ±Qo- In- 
tegrating p(k) we obtain the counting function Z(k) and 
the ground-state quasimomenta fc[/j]'s through the in- 
verse as kj = Z~ 1 (I-j/t). With kj's put in Eq. (fl~5l we 
can evaluate AEl at an arbitrary strength of interaction. 

Finally we discuss the yrast state in the light of the 
CFT. The low- lying excitation spectrum E ^ ' , among 
such eigenstates with the same I p and Ih is given by 
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Here AD represents the umklapp process with AD /I = 
(Zg(Q + ) + Ze(Q~))/2, positive integers n± denote 
particle-hole excitations near Q^, and AN = N — 
N. The Fermi velocity Vf is given hy vp — 
e (Q )/(2irpP' h (Qo)). We define dressed charge £(fc) by 

Z(k) = 1 - J% K(k - q)Z(q)dq/27r, and £ = £(Qo). 

When L is an integral multiple of N; i.e., when the 
system undergoes the center-of-mass rotation, the sys- 
tem is mapped to the L — state in a rotating frame by 
the Galilean transformation. In this case the low-lying 
excitation spectrum is consistent with those of the U(l) 
CFT with central charge c = 1. The connection also 
holds for various integrable systems with finite angular 
momentum. It enables us to predict various scaling prop- 
erties of correlation functions, such as scaling exponents. 
It should play a significant role in the study of spectral 
properties for rotating systems of interacting ID bosons 
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and in the connection to solitons 

In summary, we have constructed excited states with 
a fixed angular momentum and evaluated the excitation 
energies in terms of dressed energy through the finite-size 
asymptotic expansion. These energies include a finite- 
size correction of the order 0(l/£), which would be de- 



tected in the ultracold atomic gases. We have found 
that higher yrast excited states have quasi-degenerate 
fine structures. The special classes of angular momen- 
tum states, where L is given by an integral multiple of 
N, are related to the U(l) CFT. This relevance will pro- 
vide a powerful scheme for the analysis of the spectra, 
correlations, and scattering experiments. 
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